Properties of cage rearrangements observed near the colloidal glass transition 
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We use confocal microscopy to study the motions of par- 
ticles in concentrated colloidal systems. Near the glass tran- 
sition, diffusive motion is inhibited, as particles spend time 
trapped in transient "cages" formed by neighboring particles. 
We measure the cage sizes and lifetimes, which respectively 
shrink and grow as the glass transition approaches. Cage re- 
arrangements are more prevalent in regions with lower local 
concentrations and higher disorder. Neighboring rearrang- 
ing particles typically move in parallel directions, although 
a nontrivial fraction move in anti-parallel directions, usually 
from pairs of particles with initial separations corresponding 
to the local maxima and minima of the pair correlation func- 
tion g(r), respectively. 



PACS numbers: 61.43.Fs, 64.70.Pf, 82.70.Dd, 61.20.Ne 
Many liquids undergo a glass transition when rapidly 
cooled, where their viscosity grows by orders of magni- 
tude for only modest decreases in temperature. This 
drastic increase in viscosity is unaccompanied by sig- 
nificant structural changes; instead, the dynamics slow 
dramatically. Physically, this slowing of the dynamics 
reflects the confinement of any given particle by a "cage" 
formed by its neighbors; it is the rearrangement of the 
cage which leads to the final structural relaxation, allow- 
ing the particle to diffuse through the sample (fl. The 
dynamics of cages have been studied with scattering mea- 
surements, which probe a spatial and temporal average 
of their behavior, and with computer simulations; how- 
ever, in real systems, the actual motion of the individual 
particles involved in cage dynamics and breakup are still 
poorly understood (l]-|j|. 

In this paper, we study the motion of the individual 
particles and their neighbors during cage breakup, and 
provide the first direct experimental visualization of this 
process. We use confocal microscopy to study the motion 
of colloidal particles in a dense suspension, an excellent 
model for the glass transition |lC|-12j|. The rearrange- 
ment of cages involves the cooperative motion of neigh- 
boring particles g-0], for example as shown in Fig. ||, 
where the most mobile particles have arrows indicating 
the direction of their motion. While most neighboring 
particles move in similar directions, a significant minor- 
ity move in opposite directions, resulting in local changes 
in topology. We find particles moving in parallel direc- 
tions typically have initial separations corresponding to 
local maxima of the pair correlation function g(r), while 
pairs of particles moving in anti-parallel directions typ- 



ically correspond to local minima of g(r). We also find 
that the more mobile particles are located in regions with 
a lower local volume fraction, and higher disorder. These 
measurements provide a direct, quantitative physical pic- 
ture of the nature of cage rearrangements. 
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FIG. 1. Cut through a three-dimensional sample, with ar- 
rows indicating the direction of motion for particles with dis- 
placements Ar > 0.2 fjxa, using At* — 600 s. The sample 
has <j> — 0.52, and the cut is 2.5 /im thick (~ 1 layer of parti- 
cles) . The arrows are all the same length in three dimensions, 
so shortened arrows indicate motion in or out of the picture. 
Lighter colors indicate particles with larger displacements. In- 
set: 120 min trajectory of one particle from this sample. 

We use sterically stabilized poly- (methylmcthacrylate) 
particles with a radius a — 1.18 /an p, [l0 Jl3 1 , immersed 
in a mixture of decalin and cycloheptylbromide. This 
solvent simultaneously matches the particle index of re- 
fraction and density to mitigate the effects of scatter- 
ing and sedimentation. Hard sphere particles undergo a 
freezing transition at a volume fraction <fif = 0.494, a 
melting transition at <j) m = 0.545, and have a glass tran- 



sition at 4> g ~ 0.58 |10| , |12| . To visualize the particles, we 
stain them with fluorescent rhodamine dye; this imparts 
a slight charge to the particles, shifting the phase transi- 
tion boundaries to <j)f « 0.38, <j) m « 0.42, and (j> g sw 0.58. 
We image a volume 60 /im x 60 /im x 10 pan, containing 
several thousand particles, and identify particle centers 
with an accuracy of at least 0.05 /im |T^,|lJ]. A typical 
particle trajectory is shown in the inset of Fig. pi; it ex- 
hibits caged motion, with a sudden cage rearrangement 
which lasts ^600 s. 

The effect of cages on the ensemble dynamics is well 
known |l|-|9|,[l2[ and is evident from the particle mean 
square displacement, (Aa; 2 ), shown in Fig. 0(a) for three 
supercooled colloidal fluids (<p < <f> g ). At the earliest 
time scales, particle motion is diffusive, as they have not 
moved far enough to encounter the cage formed by their 



neighbors. As their displacement becomes larger, their 
motion is impeded by the cage, leading to the plateau 
in (Ax 2 ). The displacement at the plateau decreases 
as 4> increases, reflecting the smaller cage size. More- 
over, the cages become more long-lived with increasing 
4>; this presumably results from the fact that cage rear- 
rangements involve a larger number of particles as the 
glass transition is approached |§-0 . The cage rearrange- 
ment itself leads to an upturn in (Ax 2 ) at the end of the 
plateau on the log-log plot; at even longer lag times, the 
asymptotic motion again becomes diffusive, albeit with a 
greatly decreased diffusion coefficient, D OQl as indicated 
by the dashed lines in Fig. 0(a). 
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FIG. 2. (a) Mean square displacements. The diagonal 
straight lines indicate estimates for (Aa; 2 ) ~ 2D oc At. Ver- 
tical bars indicate At*. (Ax 2 ) is shown rather than (Ar 2 ) 
as the poorer z-resolution of the microscope artificially in- 
creases (Ar 2 ). (b) Nongaussian parameter «2(Ai). (c) The 
cage correlation function. 

We estimate the time scale for cage rearrangement, 
At* , by finding the maximum of the nongaussian param- 
eter a 2 (At) = (Ax 4 ) / {3(Ax 2 )) - 1, shown in Fig. |(b). 
a2 is computed from the ID displacement distribution 
function P(Ax; At) |^|, and is zero for a gaussian dis- 
tribution, and largest when there are broad tails. To 
compare these distributions for different 0, we normal- 
ize the displacements Ar by the average, f = (Ar 2 ) 1 / 2 . 
The radial step size distribution P(Ar/f; At*) is plotted 
for three values of in Fig. ||(a). There are more large 
displacements than expected for a gaussian distribution, 



shown by the dashed line; moreover, as <p g is approached, 
the tails become even broader, as shown by the circles for 
<f> = 0.56. This reflects the anomalously large motion of 
the particles undergoing cage rearrangements. 
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FIG. 3. (a) Step-size distribution. The vertical dashed 
line indicates the 8%ile cutoff. The dotted line is a gaussian 
with width f — 1. (b) C cag c(Ai*; Ar). (c) Average subsequent 
displacement (a;') along the direction of the original displace- 
ment Ar. The dashed line is a linear fit to the small Ar region 
for the (j> = 0.56 data (circles), (d) The average number of 
ordered neighbors N particles have. The symbols are: trian- 
gles (j) — 0.46, At* = 300 s, r = 0.61 /xm; squares <j> = 0.52, 
Ai* = 600 s, f = 0.22 ^m; circles </> = 0.56, At* = 1000 s, 
r = 0.16 nm. 

These large displacements contribute to the increase 
in the mean square displacement at lag times greater 
than At*. To confirm that this motion directly reflects 
structural relaxation, we calculate a topological cage cor- 
relation function, C ca go(At) M. We define particles as 
nearest neighbors if their separation is less than a cut- 
off distance set by the first minimum of g(r); C cagc (Ai) 
is the fraction of particles with the same neighbors at 
time t and time t + At, averaged over all t. As shown in 
Fig. 0(c), for the lowest volume fraction, <fi = 0.46 (tri- 
angles) , the particles are barely caged, but as <j> increases 
toward <j) g , the decay of C cagc slows dramatically. To 
determine which particles are responsible for topological 
changes, we replot C cago (At*) as a function of the nor- 



malized displacement. As shown in Fig. ||(b), C cago (At*) 
decreases significantly for particles with larger displace- 
ments, confirming that these particles contribute most to 
the structural relaxation. 

To directly measure the size of the cage, we investi- 
gate the temporal correlations of the motion of individual 
particles |J. Particles that remain caged must have no 
significant net displacement over long times; by contrast 
those whose cages rearrange do have net displacements. 
To quantify this, we compare a particle's displacements 
Af and Af during sequential time intervals of At*. Wc 
focus on x' , the component of Af along the direction of 
the original displacement Af; (a;') is plotted in Fig. ||(c). 
(x') is always less than zero, indicating that the average 
motion is anticorrelated. For small initial displacements 
the behavior is linear, x' — — c|(Af)|; larger values of 
c indicate more highly anticorrelated motion J9[. The 
cage constrains a particle so that the farther a particle 
moves, the farther it will move back. The linear relation- 
ship fails for larger displacements, as seen in Fig. ||(c); 
thus the degree of anticorrelation actually decreases - 
particles move shorter distances (x'} than expected by 
linear extrapolation from the small Af behavior. Wc 
identify the end of the linear regime as the cage size, 
''cage Kj. For comparison, we estimate the size of the 
cage by other methods. The simplest is based on the free 
volume of the system compared to random close packing, 
''free = 4a[((/> T . cp /(^) 1 / 3 — 1]. Another estimate comes from 
the mean square displacement at the cage rearrangement 
time scale At*, r msd = (Ar 2 (At*))5. The final estimate 
is Ar*(At*), chosen so that 5% of particles have displace- 
ments Ar > Ar* [g|| . The values of these estimates are 
listed in Table 1 and are in good agreement. The cage 
size decreases as the glass transition is approached, in 
agreement with expectations [p|,|i"5||. 

We can use this cage size to model the particle mo- 
tion as a random walk, alternating between steps (cage 
rearrangements) and pauses (caging). At long times, the 
diffusive motion reflects particles undergoing many steps 
of size r cage in random directions. ^From the Central 
Limit Theorem, fl M = r 2 ag0 /2At** , where At** is the 
mean time between steps, or the average cage lifetime 
|l6| . The measured values of D^ , and the values calcu- 
lated for At** are listed in Table 1. We find that At** is 
significantly larger than At* , the time scale for a particle 
to move during one cage rearrangement. Both grow as 4> 
increases, consistent with previous experiments |T|, pT] . [l2| 
and simulations [B-Kl|L5L|[ ■ While the cage size decreases, 
it remains finite as <j) g is approached, implying that the 
dramatic decrease in Doo is due primarily to the increas- 
ing cage lifetime. Our random-walk picture also yields 
an unambiguous estimate of the fraction of particles in- 
volved in cage rearrangements at any given time; it is 
given by the ratio At* /At**. This ratio is -8% for all 
samples except the most dilute (for which it is 16%); it 
is shown by the vertical dashed line in Fig. [|(a). 

Several computer simulations looked for the underly- 
ing origins of cage rearrangements, and found correla- 



tions between mobile particles and their local environ- 
ments P,[^ [l7lJlgll . Mobile particles tended to be in re- 
gions with higher disorder [|7| Jl7|Jl^ | , lower density [f7|,^8[, 
and higher potential energy (and thus higher disorder) 
B . Ultimately, our colloidal samples will crystallize, and 
it is possible that the evolution from disorder to crys- 
talline order drives the structural relaxation in our ex- 
periments. To investigate this possibility, we quantify 
the local order with an order parameter that identifies lo- 
cal crystalline regions from particle positions |L!| . If two 
adjacent particles have similar orientations of their neigh- 
bors about each of them, the two particles are ordered 
neighbors p0| |; most particles have between 2 and 4 or- 
dered neighbors. Particles with larger displacements have 
fewer ordered neighbors, as shown in Fig. |3j(d), in agree- 
ment with the correlation between mobility and disorder 
seen in simulations |g,0,|l7|,[18|] . Furthermore, the mobile 
colloidal particles typically move to positions where their 
number of ordered neighbors has increased by ~1, sug- 
gesting a slow evolution toward crystalline structure. We 
emphasize, however, that crystallization occurs on time 
scales significantly longer than these observations; more- 
over, crystallization is a nucleation and growth process, 
with full crystalline order of the sample resulting from 
the growth of only a small number of nuclei. Thus, any 
effect of the local crystal order on the cage rearrange- 
ment may be driven more by local variations in volume 
fraction, rather than an evolution to the state with the 
lowest global energy minimum fll9[ |. In fact, we find that 
the volumes of the Voronoi polyhedra associated with 
the mobile particles are larger on average, giving a local 
volume fraction as much as 5<f> — 0(0.03) lower than the 
mean. This suggests that the particles with smaller (f> are 
farther from <p g and thus are more likely to rearrange, in 



agreement with the simulations 

Cage rearrangement is not strictly a localized event; in- 
stead, many neighboring particles are typically involved; 
they often move in similar directions, as shown by the 
arrows in Fig. |lj |2|-|7[ . To quantify the propensity for mo- 
tion of neighboring particles in the same direction, we cal- 
culate the distribution of angles, P(9), between the dis- 
placement vectors of all neighboring particles, measured 
at At* to define the displacements. The probability of ob- 
serving two displacement vectors forming an angle in the 
range (9, 9 + d9; 4>,4> + d<p) is given by P(0, <j>) sin 9 d0 d<fi; 
in Fig. 0(a) we plot P(0, <j>) which, by symmetry, depends 
only on 9. This function is strongly peaked at « 0, in- 
dicating that two particles usually move in parallel direc- 
tions ("strings") [£|-0, although a significant fraction of 
particles move in anti-parallel directions (9 « n, "mixing 
regions") |2l]] . When 9 w tt, particles are either converg- 
ing or diverging; examples of each case can be seen in 
Fig. 0. We find that these "mixing" pairs of particles 
(9 > 7f/2) are approximately twice as likely to result in 
the two particles no longer being neighbors (compared 
to the pairs with with 9 < ir/2). Thus, these rearrange- 
ments, while less frequent, are responsible for much of 
the topological changes. 
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FIG. 4. (a) Probability distribution of angle between dis- 
placement vectors for two neighboring particles, as a function 
of the polar angle 9. (b) Pair correlation function g(Ri2) 
(solid line) and width obxi of probability distribution func- 
tion P(ARi2\Ri2) (dashed line). Insert shows P(AR 12 \Ri2) 
for nearest neighbors, with initial separations of 2.4 /an 
< R12 < 2.6 (im (solid line), and anti-neighbors with 3.0 
fim < R12 < 3.2 /xm (dotted line). All data are for <j> — 0.56 
and At* = 1000s. 

To investigate the nature of the correlated motion of 
neighboring particles, we calculate the probability dis- 
tribution function, P{AR\2\R\2), which measures the 
change in separation, Ai?i2, between two nearby parti- 
cles, initially separated by R\2, after a time interval At*. 
These distributions are narrower if R\2 corresponds to a 
peak of the pair correlation function, g(R\2), than if R\2 
corresponds to a minimum, as shown by the solid and 
dashed lines in the inset of Fig. 0(b), respectively. This is 
shown more dramatically in Fig. [|(b), where we compare 
the widths of the distribution functions, <Jr 12 , directly 
with g(Ri2), showing that they are anticorrelated. This 
is reminiscent of the behavior of the collective diffusion 
coefficient, which varies as the inverse of the static struc- 
ture factor, D(q) ~ l/S(q) |22| ]. The relaxation of fluctu- 
ations at the peak of the structure factor, which reflects 
the most favorable structure, is slowed relative to other 
values of q. Our results show a similar behavior occurs 
in real space; particles whose separation corresponds to 
a peak in g{R\2) are in more favorable relative positions 
and tend to move together, so that their separation does 
not change; by contrast, particles whose separation corre- 
sponds to minima of g{R\2) are in less favorable relative 
positions, and tend to move in antiparallel directions. 

This work reveals a physical picture of cage-trapping 
and rearrangement. Cage rearrangements involve local- 
ized clusters of particles with large displacements, in re- 
gions with higher disorder and higher free volume. Re- 
arranging particles typically move in parallel directions 
("strings"). However, there are also "mixing regions" 
where particles move in other directions [Fig. ||(a)], ac- 
counting for a significant fraction of the topological rear- 
rangements. Both the collective nature of the relaxation, 
and the local origin of the cage rearrangements, clearly 
play a key role in the behavior of supercooled fluids of 



colloidal suspensions; it remains to be seen whether they 
also are important effects in other glasses. 
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where B(i) is the number of nearest neighbors of 
particle i, C is a normalization factor such that 
J2 m qirn.(i)qfm(i) = 1> ?ij is the unit vector pointing from 
particle i to its neighbor j, and Yj m is a spherical har- 
monic function. Neighboring particles i and j are ordered 
neighbors if | }] qi m (i)q* m (.i)\ exceeds a threshold value 
of 0.5. As in other studies |l9| , we used I = 6. 
The distribution P(8, <j>) remains qualitatively unchanged 
even considering only pairs of particles with large dis- 
placements, which are involved in cage rearrangements. 
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TABLE I. Estimates for the cage size, the cage rearrange- 
ment time scale At*, the cage lifetime At**, and the asymp- 
totic diffusion coefficient Doo. 
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